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A test particle of mass (iona bound geodesic of a Kerr black hole of mass M 3> A* w iU slowly in- 
spiral as gravitational radiation extracts energy and angular momentum from its orbit. This inspiral 
can be considered adiabatic when the orbital period is much shorter than the timescale on which en- 
ergy is radiated, and quasicircular when the radial velocity is much less than the azimuthal velocity. 
Although the inspiral always remains adiabatic provided fi <C M, the quasicircular approximation 
breaks down as the particle approaches the innermost stable circular orbit (ISCO). In this paper, 
we relax the quasicircular approximation and solve the radial equation of motion explicitly near 
the ISCO. We use the requirement that the test particle's 4-velocity remains properly normalized 
to calculate a new contribution to the difference between its energy and angular momentum. This 
difference determines how a black hole's spin changes following a test-particle merger, and can be 
extrapolated to help predict the mass and spin of the final black hole produced in finite-mass-ratio 
black-hole mergers. Our new contribution is particularly important for nearly maximally spinning 
black holes, as it can affect whether a merger produces a naked singularity. 



I. INTRODUCTION 



sor Q^v [8]: 



Supermassive black holes (SBHs) with masses 
1O 6 M < M < 10 10 Af Q reside at the centers of most 
large galaxies [T]. These SBHs will be surrounded by 
dense cusps of stars [2], some fraction of which will con- 
sist of compact objects (white dwarfs, neutron stars, 
or stellar-mass black holes) of mass /i ~ 1 — 10 M@. 
Compact objects whose orbital velocities lie within a 
"loss cone" about the radial direction [3] will inspiral 
into the SBHs under the influence of gravitational ra- 
diation [3 J. These extreme-mass-ratio inspirals (EM- 
RIs) will be an important source of gravitational waves 
(GWs) for the proposed space-based GW detector LISA 
[5]. LISA will rely on matched filtering to detect EM- 
RIs, requiring templates that accurately track the phase 
of the GWs when their fundamental frequency / is in 
the range 10 -4 < / < 10 _1 Hz. These GW frequencies 
correspond to orbital frequencies f2 = nf at the inner- 
most stable circular orbits (ISCOs) of SBHs with masses 
1O 5 M < M < 1O 7 M . To maximize the number of EM- 
RIs detected by LISA, it is therefore essential to under- 
stand the trajectories of compact objects as they inspiral 
all the way to the ISCO. 

In the absence of gravitational radiation, test parti- 
cles travel on geodesies of the Kerr metric [S] that 
describes the spacetimc of spinning black holes. Kerr 
geodesies are characterized by four constants of mo- 
tion: the particle's rest mass /i, energy E, z component 
of angular momentum L z , and Carter constant Q [7j. 
These constants can be determined from the particle's 
4-momcntum and the Kerr mctric's timelikc Killing 
vector Tp, azimuthal Killing vector and Killing ten- 



* Electronic address: mhklO@nyu.edu 



/' 

E 

Q 



(la) 
(lb) 
(1c) 
(Id) 



Equatorial geodesies have Q = 0; for computational sim- 
plicity we will restrict our attention to equatorial orbits 
for the remainder of this paper. 

In the test-particle limit r\ = [i/M <C 1, the stress- 
energy tensor of a test particle moving on a Kerr geodesic 
sources gravitational radiation that can be calculated us- 
ing black- hole perturbation theory [S]. This radiation 
extracts energy and angular momentum from the orbit, 
causing the particle to migrate through the phase space 
{E,L Z }. The GW energy flux Eqw and angular momen- 
tum flux LzQw (an overdot symbolizes a derivative with 
respect to Boyer-Lindquist coordinate time) are propor- 
tional to rj 2 , implying that the timescale £gw ~ E/Eqw 
is proportional to r^ 1 . The orbital period t or b ~ is 
independent of rj to lowest order, implying that the inspi- 
ral will always be adiabatic (tew ^orb) for sufficiently 
small 7/. Gravitational radiation circularizes eccentric or- 
bits in the post-Newtonian regime |10) , implying that the 
early inspiral will also be quasicircular (r <C rfi). Dur- 
ing this adiabatic, quasicircular portion of the inspiral, 
the trajectory of the test particle is well described by a 
sequence of circular geodesies of ever decreasing Boyer- 
Lindquist radius r. In the quasicircular approximation, 
the radial velocity is given by 



E, 



GW 
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where dE/dr is the derivative of the energy of a circular 
equatorial geodesic with respect to its Boyer-Lindquist 
radius. This radial velocity is proportional to rj, and 
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approaches negative infinity as r approaches nsco where 
dE/dr vanishes. This behavior is unphysical, and reflects 
the breakdown of the quasicircular approximation as the 
test particle approaches the ISCO. 

In the vicinity of the ISCO, the test particle transi- 
tions from the quasicircular inspiral described above to 
a "captured" plunge that crosses the event horizon [TT] . 
This transition regime has been investigated in the test- 
particle limit by Ori and Thorne [12] (hereafter OT) and 
for nonspinning black holes or arbitrary mass ratio by 
Buonanno and Damour |13j . Sundararajan |14) general- 
ized the approach of OT to inclined and eccentric orbits, 
deriving trajectories that served as sources for the gravi- 
tational waveforms calculated in Sundararajan, Khanna, 
and Hughes [IS]. We will adopt the notation of OT 
throughout this paper to facilitate comparisons between 
our results and theirs. In Sec. [IT} we will review the previ- 
ous treatment of the transition regime and the discovery 
of scaling relations that can be used to apply a universal 
dimensionless trajectory to mergers with arbitrary mass 
ratio and black-hole spin. In Sec. |III| we show that the 
test particle's 4- momentum is not properly normalized 
according to Eq. (la) in this model of the transition. We 



develop a new model that satisfies this requirement by 
allowing the particle's energy and angular momentum to 
vary independently. We then explore this new model's 
predictions for the behavior of the transition regime for 
nearly maximal black-hole spins in Sec. |IV| A summary 
of our principle findings and their implications is given 
in Sec. El 



II. PREVIOUS TREATMENT 

A. Geodesic Motion 

We begin by reviewing the motion of test particles on 
Kerr geodesies. In the equatorial plane (9 — tt/2), the 
Kerr metric for a black hole of mass M and spin a can 
be written in Boyer-Lindquist coordinates [16] as 
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while according to Eqs. (lb I and (lc), the dimensionless 



energy E = E/fi, and angular momentum L = L z j^xM 
are given by 



E 



L 



2 

1- - 
f 



•r 2 



a 2 



dt 2d dip 
df f df 
To 2 \ d(f> 
f I df 



2d dt 
f df 



(5a) 
(5b) 



Solving Eqs. (5a) and (5b) for dt/df and d<p/df, inserting 
the result into Eq. Q, then solving for (df/df) 2 yields 



where 



V(f,E,L)=l- 
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is the effective potential. 

Geodesic motion is alternatively described by the 
geodesic equations 



d 2 ^ p df a dx 13 
df 2 + a P~df~df 
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where T^g are the Christoffel symbols for the Kerr met- 
ric. When E and L are constants, the second of these 
equations (x 1 = f) is equivalent to the derivative of 
Eq. ([6} with respect to f 
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B. Near the ISCO 



ds 2 = - { 1 - — ) dt 2 



2 f) 



1 - 



r 2 d6 2 



2Ma 2 



2M a z 

— + ~2 

,, AM a 



dr 2 
dtdcj) . (3) 



Here and throughout this paper we use units in which 
Newton's constant G and the speed of light c are unity. 
We can also define a dimensionless radius f = r/M, coor- 
dinate time t = t/M, and spin d = a/M. In these coor- 



dinates, the three constants of motion given in Eqs. ( la 



through (lc) provide three equations for the evolution of 
t, f, and <j) as functions of the dimensionless proper time 
f = t/M along the particle's worldline. Eq. (la I can be 



In the vicinity of the ISCO, Eq. |6} can be Taylor ex- 
panded about the ISCO values of f, E, and L. Expanding 
in the small variables 

R = f - nsco , (10a) 
X = n- s } co (E - E lsco ) , (10b) 
f = L - Zisco , (10c) 
where fiisco is the orbital frequency at the ISCO, yields 
dR\ 2 2a - dV . 

df) =- T I * + w R e + 8L<x-e> 



~ d 2 V 



(11) 
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where the ellipsis denotes higher-order terms. Following 
the OT notation we define 



= I f d3V \ 

1 / d 2 V - d 2 V \ 

p = --[ — = — + n^= — 

2\dLdf dEdfJ 
and we have made use of the relation 
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which holds at extrema of the effective potential. Again 
assuming that E and L are constant, Eq. (Ill is equiva- 
lent to 

d ^ = -aR 2 +^- 1 -n d X( X -0- (14) 
dr 2 ^ 2 dEdr 

Terms of 0(R 3 ), 0(£ 2 ,X 2 ) an d higher order have been 
neglected in Eq. (14). 




C. Radiation Reaction 

We have so far neglected the effects of radiation re- 
action on the test particle's motion. The proper way 
to include radiation reaction would be to calculate the 
conservative and dissipative parts of the self- force [T7] . 
and add these additional terms to the right-hand side of 
Eq. p]). Although self- force calculations have progressed 
rapidly in recent years (for recent reviews see [T5I [TT?] ). 
the full self-force is not yet available for test particles on 
Kerr geodesies near the ISCO. Without this self-force, 
the effects of radiation reaction can be approximated by 
making the energy and angular momentum time depen- 
dent in Eqs. ^ and 

For particles on circular orbits, the rates at which en- 
ergy and angular momentum are radiated are related by 



dE 
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dr 
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OT assume that this relation holds throughout the tran- 
sition, implying that 



where 



X = € = T\KT 



\ dt dr 
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We will later relax this assumption in Sec. |III| Inserting 
Eq. fl6|) into Eq. Jli| yields 



d2R D 2 a - 

d,T z 

identical to Eq. (3.15) of OT. 



(18) 



FIG. 1: The dimensionless radius X as a function of dimen- 
sionless time T during the transition from adiabatic inspiral to 
plunge. The solid black curve shows the numerical solution to 
Eq. (211, while the dashed blue and dotted red curves show 



the approximate analytic solutions at early and late times 
given by Eqs. (22 1 and (231 respectively. This figure is a re- 



production of Fig. 2 in OT. 



D. Dimensionless Equation of Motion 



Previous studies [T2J[T3] noticed that Eq. (18) can be 
converted into dimensionless form by defining 



R = r) 2/5 R X , 



where 
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In these variables, Eq. ( 18 ) becomes 
d 2 X 



dT 2 



T . 
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At early times, the particle's radial velocity and acceler- 
ation approach zero, suggesting that the correct solution 
to Eq. (21) asymptotes to 



T 



(22) 



At late times, the particle plunges into the horizon (X = 
— oo ) in a finite proper time f. The second term on the 
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right-hand of Eq. ( 21 ) can then be neglected, yielding the 



approximate solution 
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Numerically, T p i ungc = 3.412. The numerical solution to 



Eq. (21 1 and approximate analytic solutions of Eqs. (22) 



and 



and (23 1 are given by the solid black, dashed blue, 
dotted red curves in Fig. [I] 

In addition to the approximate radial trajectory X(T), 
this model also provides an estimate of the energy and an- 
gular momentum radiated during the transition. Assum- 
ing that energy and angular momentum are radiated at 
the ISCO rate throughout the transition, Eq. (23) implies 



that the test particle will radiate an additional amount 



AE tI = Eisco 
AL tl - = Ligco 



£Wai = ?7 4/5 ^iscoKToTp lungo (24a) 

£final = 77 4/5 KT T plungc (24b) 



beyond that calculated in the adiabatic approximation. 
Since AE tI and AL tr scale as rj 4 / 5 , in the test-particle 
limit rj — > the specific energy E and angular momen- 
tum L radiated during the transition should exceed that 
radiated during the ringdown which scales as rj. Given 
the success of efforts to predict black-hole spins by ex- 
trapolating from the test-particle limit f2"0H2"3"] , a closer 
examination of the energy and angular momentum radi- 
ated during the transition is worthwhile. 



III. NORMALIZATION OF 4- VELOCITY 



We noted in Sec. Ill Al that the normalization of the 4- 
velocity Q and the geodesic equation ([8| provide alter- 
native equations of radial motion. While they are equiv- 
alent for geodesies, they differ once E and L become time 
dependent. Ori and Thorne [12] solve the geodesic equa- 
tion under the assumption that E and L vary with proper 
time according to Eq. (16). What does this assumption 



imply for the norm of the 4-vclocity? Expressed in the 
dimensionless variables of Eq. (19), Eq. (11) becomes 



dX_ 



= --X'- 



2XT . 
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The left and right-hand sides of Eq. ( 25 ) for the numerical 
solution X(T) of Eq. ^ are plotted in Fig. [5] The 
solid black curve, showing the left-hand side (dX/dT) 2 , 
is positive definite as one would expect. The dashed blue 
curve shows the right-hand side. If Eq. ( 25 1 was satisfied, 



the two curves would be identical and the right-hand side 
would be positive definite as well. This is clearly not the 
case. 

What docs this mean physically? According to 
Eq. (16), the dimensionless variable T is simultaneously 
proportional to the proper time f, the energy x, an d 
the angular momentum £. Values of X(T) for which 
the right-hand side of Eq. (25) is negative correspond 




FIG. 2: The square of the dimensionless radial velocity 
dX/ dT as a function of dimensionless time T during the tran- 
sition from adiabatic inspiral to plunge. The solid black curve 
shows the left-hand side of Eq. ( 25 I for the numerical solution 



to Eq. ( 21 1, whi le the dashed blue curve shows the right-hand 
side of Eq. ( 25 1 for this same numerical solution. The two 



curves clearly differ, and the dashed blue curve becomes neg- 
ative when it falls below the dotted red line. Since (dX/dT) 2 
is positive definite, this indicates that the numerical solution 
occupies a forbidden portion of {X, T} parameter space. 



to values of {f, E, L} for which the effective potential 
V(f, E, L) exceeds the square of the energy E 2 in viola- 
tion of Eq. ([6]). For example, at T — the energy and an- 
gular momentum equal their ISCO values, yet X(0) > 
as seen in Fig.[j] Since V(fisco, -E"iSCO, -^isco) = Ef sco 
and the effective potential V is a monotonically increas- 
ing function of X for T — 0, the right-hand side of 
Eq. ( 25 ) must be negative at T — as seen in Fig. [2] 
The discrepancy between the left and right-hand sides 
of Eq. ( p5| ) is not a consequence of a breakdown in the 
Taylor expansion of the effective potential for R, £ > 1. 
Figure[2]shows that Eq. ( 25 ) is violated even for X, T ~ 1, 
which corresponds to R, £ <C 1 when ij<1. 

How can we reconcile the normalization of the 4- 
velocity Q and the geodesic equation ( f8|) i n the presence 
of radiation reaction? Although Eqs. (|15|) and (16) cor- 



rectly specify the energy, angular momentum, and their 
derivatives on a circular orbit at ngco; there is no reason 
to expect them to hold exactly when the radial velocity 
is nonzero as it is during the transition. If we relax the 
requirement that x = d we recover the third and fourth 
terms on the right-hand side of Eq. ( 11 ) that vanished in 



the previous treatment. We can retain the dimensionless 
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FIG. 3: The solid black, short-dashed blue, and dotted red 
curves in this plot are identical to those in Fig. [2] but we 
have also added the dimensionless difference between the en- 
ergy and angular momentum Y(T) which is given by the long- 
dashed green curve. The short-dashed blue c urve shows the 
first two terms on the right-hand side of Eq. (28 1, while the 



long-dashed green curve shows our new third term. Their sum 
equals the solid black curve, implying that Eq. | ]28| is now sat- 
isfied and the test particle's 4- velocity is properly normalized. 



form of the equations of motion by defining 



where 



(X-Oo 



~ 4 / 5 (/3«) 6 / 5 



dV 
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(27) 



The second-order equation (21) remains unchanged to 
lowest order in 77, while Eq. (25) gains an additional term 
to become 



dX 



-X 3 



2XT + Y 



(28) 



Instead of Eq. ( 15 ), which would imply that Y is constant 



throughout the transition, we evolve Y according to 

2X . (29) 



dY 



This is precisely what is required to restore the consis- 
tency between Eqs. (21 ) and (28) that exists for geodesic 



motion. Since Eq. (21) remains unchanged, its solution 



X{T) remains unchanged as well. 



We can solve Eq. ( 29 ) by inserting X(T) into the right- 



hand side and choosing the correct initial condition. This 
initial condition can be found by matching to the qua- 
sicircular inspiral at early times. During the inspiral, E 
and L are given by [TT] 



1-2/f + a/f 3 / 2 
v/l-3/f + 2a/r 3 /2 ' 
= , 1/2 l-2a/f3/3 + s2/r - 2 
v/l-3/r + 25/r 3 / 2 
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L, 



Taylor expanding about the ISCO, 

X~t = n^ co (E c (r) - E lsco ) - (L c (r) 

d 3 L r \ 
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Using the definitions of X and Y in Eqs. ( 19a) and (26 1, 
this implies 



X 3 



6 ^ df 3 ^ 3 y isco Ui/ IS co 
= -\x 3 T — > —00 . (32) 



(28 



Inserting Eqs. \22\ and (|32|) into the right-hand side of 
Eq. 



shows that dX/dT vanishes as T — > — 00 as 
required by Eq. ( [22] ) . The numerical solution Y(T) with 
this initial condition is shown by the long-dashed green 
curve in Fig. [3] With the addition of this new term to 
the right-hand side, Eq. (28) is now satisfied ensuring 



that the 4-velocity is normalized according to Eq. ( la 



Our new term Y(T) is important for more than just 
mathematical consistency. Equation (26) suggests that 
the physical quantity x — £ is proportional to rj 6 ^ 5 dur- 
ing the transition, making it higher order than self-force 
corrections that scale as 77. However, as the test particle 
plunges into the horizon, the asymptotic solution (23) 
and Eq. ( 29 ) imply 



Y(T) 



f dY 



Y(0)- 



12 



(-^plunge E 



■l\2 
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12 



plunge 



(33) 



which diverges as T — ¥ T p i ung0 . This divergence can be 
seen from the limiting behavior of the three curves in 
Fig. |3j Although the solid black and short-dashed blue 
curves appear to converge as T — > T p i ungc , their difference 
Y(T) shown by the long-dashed green curve in fact di- 



verges in accordance with Eq. ( 33 1 . Does the divergence 
of the dimensionless Y(T) impl 
the physical quantity x. 
into the black hole? 



dy a similar divergence in 
£ as the test particle plunges 
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To answer this question, we must examine the validity 
of the asymptotic solution (23) as T — >• T p i unge . Footnote 



3 of OT notes that the divergence of X(T) results from a 
breakdown in the dimensionless equation of motion (21) 



when higher-order terms in the Taylor expansion about 
the ISCO become important. This breakdown occurs at a 
coordinate radius i?brcak> which corresponds to a dimen- 
sionless radius Xb rC ak oc r]~ 2 / b according to Eq. (19a). 
This in turn implies that T p i ungG — Tbrcak oc f] 1 ^ by 
Eq. (23) and Ybrcak °c rj^ 1 ^ from Eq. (33). According to 



Eq. (26), the test particle's energy during the transition 



will receive a correction 



AE r: 



V 6/5 ^isco(x-OoY 



(34) 



When Y ~ Ybrcak, this correction will be linearly pro- 
portional to 1) just like self-force corrections. It will be 
instructive to compare future self-force and time-domain 
perturbation theory calculations with this analytic result. 



IV. MAXIMAL SPINS 

Kerr black holes [6] have spins 5 < 1; objects with 
larger spins are "naked singularities" unclothed by an 
event horizon. Penrose [21] proposed that a cosmic censor 
protects general relativity by preventing the formation of 
such naked singularities, but this cosmic censorship con- 
jecture has never been proven in full generality. Wald [55] 
determined that maximally spinning black holes could 
not accrete test particles that would drive them over the 
Kerr limit 5 = 1, but Jacobson and Sotiriou |26j recently 
showed that black holes with spins 8 = 1 — a <C 1 could 
be spun above this limit by accreting particles with a 
finite energy E = r/ME. Barausse et al. [27] showed 
that gravitational radiation could not prevent some of 
these finite-mass-ratio mergers from producing a naked 
singularity, but that self-force corrections could be signif- 
icant. Given that finite-mass-ratio effects may determine 
whether naked singularities can exist, it is worthwhile 
to examine what happens to the transition region in the 
maximally spinning limit. 

The previous treatment of the transition from quasi- 
circular inspiral to plunge assumes that the effective po- 
tential V and fluxes dE/df and dL/df can be Taylor 
expanded in a neighborhood of the ISCO. However, as 
the black hole's spin approaches the Kerr limit a — > M, 
r iSCO an d the outer horizon r + both approach M in 
Boyer-Lindquist coordinates. This suggests that it may 
not be possible to construct a neighborhood of the ISCO 
that does not include the horizon. To determine whether 
or not such a neighborhood exists, we must compare 
R + = hsco — f+ with i? in Eq. (20a I as 5 — > 0. In 
this limit ITTl. 
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FIG. 4: The general-relativistic correction £ to the Newto- 
nian quadrupole-moment formula for the GW luminosity as 
a function of black-hole spin 5 = 1 — 5. The points have 
been calculated using the black-hole perturbation theory code 
GREMLIN [30] and fitted with a power law £ = A8 m . 



The behavior of Rq as S — > depends on the behavior 
of the fluxes dE/df and dL/df. These in turn depend 
on £ , defined by 



•^GW — 



32 



(36) 



as the general-relativistic correction to the Newtonian 
quadrupole-moment formula for the GW luminosity [25] . 
The ISCO values of this correction £ were calculated for 
spins —0.99 < a < 0.999 in [2H], and have been calculated 
down to 5 = 10~ 4 [35] using the GREMLIN (Gravita- 
tional Radiation in the Extreme Mass Ratio Limit) code 
presented in [30]. We have fitted these calculated values 
to a power law £ = AS rn as shown in Fig. |4j The best-fit 
parameters for this power law are A = 1.80, m = 0.317. 
The summation over spheroidal harmonics of the Weyl 
scalar ip4 needed to compute £ converges very slowly in 
the limit 8 — > 0, so our best-fit parameters should be 
regarded with caution until they can be confirmed by a 
technique better suited to this limit. Chrzanowski [31] 
estimated that 



E GW ~ V 2 R 



+ ■ 



(37) 



implying R + cx 8 1 / 3, . 



which also suggests that m ~ 1/3. The rough agreement 
between our numerical fit and Chrzanowski's estimate 
gives us some confidence that m ~ 1/3 is close to the 
correct value. 
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The scale Rq of the transition region depends on the 
relativistic correction £ through k, which was defined in 
Eq. |l7|) and can be expressed as 



5? (n 7 / 3 —e 

5 V df 



ISCO 



For circular equatorial Kerr geodesies, as 5 — > 0, 



di 
df 



(45) 



-1/3 
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dV 
dL 



1 

2 

< 

16 

1 , 

4 



:j45 m-l/3 



(4<5) 



1/3 



implying 



R oc (5 2m/5 , 
t oc (T m/5 , 
(X-Oo oc 5<W5-V3 



(38) 

(39a) 
(39b) 

(39c) 
(39d) 
(39e) 
(39f) 



(40a) 
(40b) 
(40c) 



from the definitions in Eqs. (|20| and (|27|) 

t: 



According to Eq. (19a 



le transition region will be 



cut off by the horizon at a dimensionless radius 



X, 



R 



2/5fl± ^ _ ?? -2/5 ( jl/3-2m/5 ^ 

R 



as S —> 0. We plot X+(<$) for several mass ratios and our 
best-fit value of m in Fig. [5] Since m < 5/6, X + — > 
in the limit of large spins. One should not expect GW 
emission beyond X + , so Eq. ( [24] ) overestimates the ad- 
ditional energy and angular momentum radiated during 
the transition. In the limit of large spins, A_E tr and Ai tr 
will be suppressed by a factor To/T p i ungc ~ 0.21, where 
X(T o )~0. _ 

Despite this suppression, the extra energy and angular 
momentum radiated during the transition remain propor- 
tional to kt o cx (5 4m / 5 - 1 / 3 j which diverges for m < 5/12. 
Our numerical fit shown in Fig.[4]suggests that m satisfies 
this inequality. However, this divergence might not be 
physical; it could merely reflect a breakdown of the Tay- 



lor expansion of the equation of motion given in Eq. ( 14 ) 



The first two terms on the right-hand si de o f this equa- 
tion are proportional to rj 4 ^ 5 , while Eq. (26) shows that 
the third term is proportional to r/ 6 ' 5 . This higher-order 
dependence on r\ justifies neglecting this third term for 
modest spins, but we must also consider how each term 
depends on S in the maximally spinning limit. In this 
limit, the first two terms are proportional to <5 4m / 5 ac- 



cording to Eqs. (39) and (40) while the third term is 



X 




-70 



0.0001 
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FIG. 5: The dimensionless horizon radius X+ = 
—rj~ 2 ' 5 R+/Ro as a function of S = 1 — 5. The dotted red, 
dashed blue, and solid black curves correspond to mass ratios 
f] = 1(T 4 , 1(T 5 , and 10~ 6 respectively. X+ -> as S 
because the relativistic correction £ to the energy radiated in 
GWs has a power-law index m < 5/6. 



proportional to J 6 " 1 / 5 - 1 / 3 . For m < 5/6, this third term 
will dominate over the first two terms, indicating a break- 
down of the original dimensionless equation of motion 
(21). For completeness, we note that we have neglected 



an additional term 



1 fd 4 V 



12 V df 4 



R A 



(42) 



ISCO 



at C(?7 6 / 5 ) on the right-hand side Eq. (14). This term is 
proportional to <5 6m / 5 as 8 — > and thus can always be 
neglected. 



Including the third term in Eq. ( 14 ) on the right-hand 
side of Eq. (plj) yields 



d 2 X 



-X 2 -T + r/ 2/5 CY , 



(43) 



where 



C^--a- 3 / 5 (/3 K ) 2 / 5 n & T- ( 

2 yH ' dEdf V 



dV_ 



(44) 

diverges as S — > for m < 5/6. Y must evolve according 
to 



— = 2X + 2 V 2 / 5 CY — 
dT 1 dT 



(45) 




FIG. 6: Upper panel: The dimensionless radius X as a func- 
tion of dimensionless time T in the limit S = 1 — a — > 0. The 
solid black curve is the OT solution, while the long-dashed 
blue, short-dashed green, and dot-dashed red curves show our 
solutions with n = 10~ 3 and 8 = 10~ 2 , 10~ 4 and 10~ 6 respec- 
tively. Lower panel: The dimensionless difference Y between 
the energy and angular momentum for the solutions shown in 
the upper panel. 



to p rese rve the normalization the 4-velocity given by 
Eq. (28). The product n 2/5 C thus measures the devi- 



ation of X(T) from the OT solution in the limit of large 
spins. We plot our new solutions X(T, n 2 ^ 5 C) for a mass 
ratio rj = 10 -3 and spins 8 — 1CP 2 , 1CP 4 and 10 -6 in 
Fig. [6j Although our solutions diverge from the OT solu- 
tion as ?y 2 / 5 C increases, the test particle crosses the ISCO 
X = (shown by the horizontal dotted line) at a later di- 
mensionless time T. This implies that more energy and 
angular momentum are radiated and the divergence in 
AE tI and AL tr as 8 — > cannot be avoided. One must 
look to finite-mass ratio effects beyond the scope of this 
paper to eliminate this unphysical result. 

One such effect that might provide a solution to this 
problem is the spinning down of the black hole due to 
the extraction of angular momentum by the superradiant 
scattering of GWs earlier in the inspiral [23] . This effect 
implies that even an initially maximally spinning black 
hole will have 8 oc n by the time the test particle reaches 
the transition region. The additional energy and angu- 
lar momentum A£ tr ,Z tr oc ^/s^m/s-i/s K ^m/5+7/16 
would then remain finite even if 8 = initially. 

How does the energy and angular momentum radiated 
during the transition affect the spin of the final black 
hole produced in the merger? If we assume that energy 



and angular momentum are conserved after the end of 
the transition, the final spin will be 



«/ = 



r,(L 



ISCO 



[1 + v(Eisco - A£ tr + A£ norm )p 
a + ?7(Aaisco + Aa tr + Aa norm ) 



(46) 



for 77 <C 1, where in the limit 8 — > 



Aaisco 



-(45) 2 / 3 , 



(47a) 



Aa tr 



Aa„ 



-isco - 2£isco ' 
-(1 - 2f2isco)AZ t r oc rf'H Am ^, (47b) 

-2A£; norm OC 77 6/5 ( j6 m /5-l/3_ ( 47c ) 



Although our new correction Aa norm to the black hole's 
final spin is subdominant in 77, for m < 5/6 it becomes 
the dominant correction as 8 — > 0. Since Ai? norm > 
when the particle crosses the horizon in the limit 8 — > 
0, our analysis suggests that gravitational radiation in 
the transition region will not promote the formation of 
naked singularities. This supports our earlier result in 
|23j . where we showed that the superradiant scattering 
of GWs emitted at f > risco would also reduce the final 
spin for a > 0.998. 



V. DISCUSSION 

The existence of an ISCO is one of the most distinc- 
tive features of a black hole. Although the presence 
of an ISCO indirectly affects the luminosity and spec- 
tra of accreting black holes, it is more cleanly probed 
by the GWs emitted as test particles (compact objects 
like white dwarfs, neutron stars, and stellar-mass black 
holes) plunge into supermassive black holes. Such GWs 
are a primary source for the proposed space-based GW 
detector LISA. The LISA detection strategy relies on con- 
volving observations with theoretically determined GW 
templates. Understanding how a test particle's position, 
energy, and angular momentum evolve near the ISCO is 
an important first step toward constructing these tem- 
plates. 

Previous studies O [13] identified a transition region 
near the ISCO where neither the quasicircular approxi- 
mation nor the assumption of geodesic motion are valid. 
They solved the radial equation of motion in this re- 
gion by Taylor expanding it about the ISCO. Ori and 
Thorne (2000) [12 fixed the energy and angular momen- 
tum fluxes to their ISCO values, while Buonanno and 
Damour (2000) [TB] worked with the Hamilton equations 
of motion and did not need to specify an explicit energy 
flux. The dimensionless solution X(T) obtained by both 
groups can be rescaled to determine the evolution of the 
test particle's position, energy, and angular momentum 
for arbitrary mass ratios and spins. The simplicity of 
their approach and the universality of their solution are 
highly appealing, but a closer examination reveals that 



9 



their solution does not properly normalize the particle's 
4-momentum. 

We undertook this study to see whether this technical 
problem could be easily remedied, or was a symptom of 
a more serious flaw in their approach. We found that by 
introducing a correction Ai? norm to the particle's energy 
at a higher order in the mass ratio rj <C 1 , we could prop- 
erly normalize the 4-momentum without altering their 
universal solution X(T) at lowest order. This relatively 
modest correction increases our confidence in this parti- 
cle trajectory and validates its use as a source for the 
construction of GW waveforms as in j!5) . 

In addition to its role in constructing GW templates, 
gravitational radiation during the transition also affects 
whether a test-particle merger can produce a naked sin- 
gularity by increasing a black hole's spin above the Kerr 
limit 2 = 1. Our calculation of the energy and angular 
momentum flux at the IS CO as 5 = 1 — 2 — > sug- 
gests that the total energy A_E tr and angular momentum 
AL tr radiated during the transition diverge in this limit. 
This divergence cannot be physical, and must therefore 
be moderated by high-spin corrections beyond the scope 
of this paper. Despite this divergence, the change A2t r in 
the black hole's spin due to this radiation remains finite 
in magnitude and negative in sign. It therefore reduces 
the likelihood that a naked singularity will be produced, 
as does our new correction Aa norm which is also negative 
and becomes dominant as S — > 0. 

Throughout this paper, we have neglected self-force 
corrections to the particle's energy even though they are 
formally lower order in r\ than our correction AE noim oc 



r; 6 / 5 . Far from the ISCO, where \Y\ 3> 1, our correc- 
tion should dominate self-force effects despite its scaling 
with rj. Near the ISCO however, self-force corrections 
should be significant. Reliable calculations of the self- 
force do not yet exist near the ISCOs of highly spinning 
black holes. Comparing such self-force calculations with 
our analysis as they become available is a subject for fu- 
ture work. Another important future test of our analysis 
would be to use time-domain perturbation theory to cal- 
culate the energy and angular momentum radiated by a 
particle whose position is given by an appropriate rescal- 
ing of the solution X(T). Agreement with our predicted 
solution Y(T) would strongly support our conjecture that 
proper normalization of the 4-momentum can be used to 
predict GW emission. Comparison with numerical rel- 
ativity may soon be possible as well, since simulations 
with mass ratios as small as 7/ = 0.01 have recently been 
performed 132] . Our predictions could be even further im- 
proved by including the gravitational radiation emitted 
near the horizon |33| . We look forward to comparing our 
calculation to these alternative theoretical approaches in 
the near future, and comparing all these predictions to 
GW observations in the hopefully not too distant future. 
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